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I. INTRODUCTION 



TMD evolution of the Sivers asymmetry 

Daniel BociQ 
University of Groningen, Zernikelaan 25, NL-9747 A A Groningen, The Netherlands 

(Dated: April 22, 2013) 

The energy scale dependence of the Sivers asymmetry in semi-inclusive deep inelastic scattering 
is studied numerically within the framework of TMD factorization that was put forward in 2011. 
The comparison to previous results in the literature shows that the treatment of next-to-leading 
logarithmic effects is important for the fall-off of the Sivers asymmetry with energy in the measurable 
regime. The TMD factorization based approach indicates that the peak of the Sivers asymmetry falls 
off with energy scale Q to good approximation as 1/Q 0,7 , somewhat faster than found previously 
based on the first TMD factorization expressions by Collins and Soper in 1981. It is found that 
the peak of the asymmetry moves rather slowly towards higher transverse momentum values as Q 
increases, which may be due to the absence of perturbative tails of the TMDs in the presented 
treatments. We conclude that the behavior of the peak of the asymmetry as a function of energy 
and transverse momentum allows for valuable tests of the TMD formalism and the considered 
approximations. To confront the TMD approach with experiment, high energy experimental data 
f- ^ from an Electron-Ion Collider is required. 
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The Sivers effect [TJ [2] is a left-right asymmetry in the transverse momentum distribution of unpolarized quarks 
inside a transversely polarized proton. It is a kx x St correlation in the quark transverse momentum (kx) 
distribution with respect to the transverse polarization (St) of the proton moving in the ^-direction. It was first 
defined as a transverse momentum dependent parton distribution (TMD) in [3J. In [5] it was shown that the 
Sivers effect in semi-inclusive DIS (SIDIS) leads to a single transverse spin asymmetry Ajjt ~ sin((f)h — (f)g)fy T Di, 

i— —i where f^~ T denotes the Sivers effect TMD and Dy the ordinary unpolarized fragmentation function. The azimuthal 

angles 4>h and <ps are of the observed outgoing hadron's transverse momentum and of the transverse spin vector 
of the proton, respectively, in frames where the proton and virtual photon are collinear and along the z-direction. 
Such an asymmetry has been clearly observed in the SIDIS process by the HERMES [5] and COMPASS [6] 

qq experiments. Since these two experiments are performed at different energies (Q), as will future experiments at 

Jefferson Lab and at a possible Electron-Ion Collider (EIC), it is important to study the evolution of the Sivers 

l/~) asymmetry with energy scale Q. 

—J The evolution is dictated by TMD factorization [7H9]. Explicit expressions to order a s allow to obtain the 

leading order scale dependence of the cross section and its asymmetries. Evolution of the Sivers TMD and its 

fT) SIDIS asymmetry has recently been studied numerically in [10 (to be referred to as AR), [TTj . and [T2] (APR). 

t-H Earlier numerical studies of the evolution of the closely related Collins effect have been done in [13] (B01) and 

[2] (B09) , the results of which can be carried over to the Sivers effect case upon trivial replacements. The main 
difference in the approaches is that AR and APR are based on the most recent TMD factorization expressions 
of [5], whereas B01 and B09 were based on the original TMD factorization expressions by Collins and Soper in 
1981 7 . Although the new form of TMD factorization is preferred on theoretical grounds, because it takes care 
of several problematic issues with earlier forms (specifically, infinite rapidity divergences and divergent Wilson- 
line self-energies that should cancel in the cross section), it is not a priori clear that numerical results based on 
earlier expressions are invalidated or affected substantially. Especially for the limited energy ranges accessible in 
experiments previous results may still be of value and in order to judge that, the example of the Sivers asymmetry 
is considered. 

The conclusion of B01 about the energy scale dependence of the Collins single spin asymmetry (SSA) in SIDIS 
and by extension of the Sivers SSA, was that in the range Q — 30 GeV to Q = 90 GeV, it falls off approximately 
as l/\/Q (actually a range 1/Q 05 -1/Q ae was obtained numerically, upon variation of the nonperturbative input). 
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This moderate fall-off has been recently contested by APR, although no direct comparison of the same quantity 
was made. The intention here is to shed light on the differences between the approaches followed in B01 and 
B09 and the recent analyses by AR and APR, by studying one particular analyzing power expression. Although 
the various approaches all coincide in the double leading logarithmic approximation, in which the running of 
the strong coupling constant a s is neglected, this approximation is valid only in a very small range of Q values. 
The running of a s will have to be included to evolve the Sivers asymmetry from for example the HERMES and 
COMPASS Q scales ((Q 2 ) = 2.4 GeV 2 and (Q 2 ) = 3.9 GeV 2 , respectively) to scales relevant for an EIC (with Q 
values up to about 100 GeV). 

One conclusion of this paper will be that the approach of AR/APR yields to good approximation 1/Q 07 for 
the fall-off of the peak of the asymmetry, when using the nonperturbative Sudakov factor of AR in the range 
from Q ~ 3 GeV to Q ~ 100 GeV, in other words is not too far from 1/Q 05-0 - 6 found in B01 and even closer to 
what one finds using the approach in B09 (section lv|, despite their considerably different expressions. We expect 
comparably moderate modifications to apply to the Collins effect asymmetries discussed in B01 and B09. 



II. THE SIVERS ASYMMETRY AT TREE LEVEL 



The expression for the single transverse spin asymmetry Ajjx ~ sin(0^ 
convolution integral: 
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where the sum runs over all quark (and anti-quark) flavors, e a denotes the quark charge in units of the positron 
charge, and the ellipses denote contributions from other, T-even TMDs [15]. The cross section is differential in 
the invariants x — Q 2 /(2 P ■ q), z = P ■ Ph/P • Q, y = {P • q)/(P • I) ~ q/l, f° r incoming hadron momentum P, 
outgoing hadron momentum P% and beam lepton momentum I and the momentum q of the virtual photon, defining 
q 2 = —Q 2 , and differential in dfl — 2dyd<j) and d 2 q T where q T = q + x P — Ph/z, such that Q\ = q\ = —q 2 
and qx has only transverse components in the frames where the hadrons are collinear. In the frames where the 
proton and the photon are collinear, the perpendicular component of Ph satisfies: P 2 ± — z 2 Q\. At tree level the 



convolution integral with weight function w is given by: 

T[w(p T ,k T )fD} = J d 2 p T d 2 k T 6 2 (p T + q T - k T )w (p T ,k T ) f a (x,p 2 T )D a (z,k 2 T ) . 
For the evolution study we will however consider the Fourier transformed expressions, i.e. 

' [/l Dl] = / Wf elbQT f ~" {x > ^ ^ (Z ' &2) = h J dhhJ ^ b ^T) fi(^ b 2 ) D?(z, b 2 



T\ 



T 



q T p T f ± n 



-2i 

1 



2ttM 



(2tt) 2 M IlT [ ' ^i^'° > 

dbb 2 MbQr) f^ a ( x ,b 2 ) DKz^ 2 ), 



(2) 



(3) 



(4) 



where /i(x, b) and D\(z, b) denote the Fourier transforms of the unpolarized TMD distribution function /i(x,p 2 ) 
and TMD fragmentation function Di(z,k T ). We have also defined 



This yields for the analyzing power of the sin(0^ — <$>§) asymmetry 
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Keeping in mind that the flavor indices in numerator and denominator are part of separate summations, we define 

A a b(x,Z,QT) = — =; =-, (7) 

MQ T JdbbJ Q {bQ T )ft(x,V)D\(z,W) 

At tree level this would be the relevant asymmetry quantity at all energy scales. Beyond tree level it would be a 
valid expression at one particular scale only. The expressions at other scales can then be obtained by evolution 
of the parameters involved. To discuss this in more detail, we will first discuss the relevant aspects of TMD 
factorization. 

III. SCALE DEPENDENCE OF THE TMD FACTORIZED CROSS SECTION 

The proof of TMD factorization of processes such as semi-inclusive DIS or Drell-Yan, has recently been finalized 
[21 [IB]. It involves a new definition of TMDs that incorporates the soft factor, which then no longer appears 
explicitly in the cross section expression. In this TMD formalism the differential cross section of for instance the 
SIDIS process at small Q T /Q 2 is written as 



da 



dxdydzd ^ T =J dH e ' ibqT ^ q->*>v>*)+° (Qt/q 2 ) ■ ( 8 ) 

The integrand for unpolarized hadrons and unpolarized quarks of flavor a is given by 

W(b,Q;x,y,z) = £ /f (x,b;{ F ,n)D<?(z,b;( D ,n)H (y,Q;n) , (9) 

a 

The partonic hard scattering part H, which for the choice \x = Q takes the form 

H(Q;a s (Q)) ex e 2 a (l + a s (Q 2 )F 1 + 0(a 2 s )) , (10) 

where F\ denotes a renormalization-scheme-dependent finite term. The dependence of the TMDs on C,p , £d and 
/i will be discussed in detail next. Once the factorization expression is given, with all its scale dependence, the 
evolution of TMD cross sections can be obtained. To obtain the evolution of spin asymmetries in these cross 
sections, one has to include spin dependent TMDs, whose Fourier transforms can be odd under b — > —b, such as 
the Sivers or Collins effect. 

Note that Eq. ([8]) contains no integrals over the partonic momentum fractions, which only appear in the large 
Qt (or equivalently, small b) limit. What is considered large Qt depends on Q. For asymptotically large Q, the 
main contribution to the b integral is from small b values, which means that the b dependence of the (Fourier 
transformed) TMDs can be calculated entirely perturbatively. For the Q values of HERMES and COMPASS, 
this is certainly not the case. The peak of the Sivers asymmetry for Q values in the range 3 to 100 GeV will be 
located at Qt values for which the 6-integration receives important contributions from b values that do not allow 
for a perturbative calculation of the b dependence. This will require separate treatment of the small and large b 
regions. 

A. Sudakov factor 

In Eq. (|9j the Fourier transformed TMDs /1 and the hard part H have a dependence on the renormalization 
scale \i. It will be chosen [i = Q, such that there are no \nQ//j, terms in the hard part, as in Eq. ( |To| ). In order to 
avoid large logarithms, the TMDs will be taken at either the scale /u& = C\jb = 2e m jb [C\ sw 1.123), or at the 
fixed scale Qq which is to be taken as the lowest scale for which perturbation theory is expected to be trustworthy 
(a common choice is Qq — 1.6 GeV). Evolving the TMDs from the scale Q to fi^ or Q will result in a separate 
factor, called the Sudakov factor, which will be discussed in this subsection. We will primarily focus on the fixed 
scale case, such that the TMDs are always considered at the same scale Qo when integrating over b. This has the 
advantage, as will become clear below, that the remaining 6-dependence of the TMDs is perturbatively calculable. 
The fixed scale option was already suggested by Collins & Soper [7] and explicitly used by Ji et al. [HI E] , who 
called it /i/,, and in B09 and APR. It is not necessarily always the optimal choice though, which depends on 
the absence of large logarithmic corrections. For the specific quantity and the energy and momentum region 



considered here, it appears to be an appropriate choice, that moreover will allow us to compare the approaches 
of APR and B09 more directly. 

The TMDs also depend on Cf(d)i which are defined as |11) : 

( F = M 2 x 2 e 2 ( yp - y °\ ( D = M 2 e 2iy ^ Vh) /z 2 , (11) 

where ypih) denotes the rapidity of the incoming proton and outgoing hadron and the dependence on the arbitrary 
rapidity cut-off y s cancels in the cross section, where only the product CfCd ~ Q 4 enters. 

The evolution of the TMD / in both £ and /i is known and given by the following Collins-Soper and Renor- 
malization Group equations, respectively [9]: 

dln/M;C,M) = R{h .^ (12) 



dlnf(x,b;(,[i) 
din fi 



7 f( ff ( M );C/M 2 ), (13) 



where dK/dln^i = —"fK{g{^)) and 7f(s(m); C/a^ 2 ) — 7f(<7(m))1) — \lK{g{n)) 1h(C/a* 2 )- With these evolution 
equations one can evolve the TMDs to the scale \i b or Q , i.e. 

/(x,6;C,Q) - f(x,b^l^ b )exp(\n(^\K(b, f i b ) + J ^ [ 7 ^ 5 ( M ); 1) - In (^ W( 5 (^))] }, (14) 



'in, 
or 



f(x,b 2 ;(,Q) = f(x,b 2 l Q 2 ,Q )exp{]n(J^jK(b,Q )+J ^[j F (g(^ 1) - In (^f) J K (g(»))] }• (15) 

The latter expression is however not optimal 1 , since it does not take care of possible large logarithms in fib/Qo- 
It is more appropriate to use instead [Hj: 



f(x,b 2 ;(,Q) - f(x,b 2 ;Qi,Q )exp\ In ( ^ 1 lK(b lf i b ) + J -^-f K (g(n)) 



[ 7F (g(^l)-ln(^j 7K (g(iJ,))]\. (16) 



d(x | 

Similar equations can be obtained for the TMD fragmentation functions D, cf. [TU]. This yields the expressions: 
^{x,b 2 ] Q F ^)b\{z,b 2 ] Q D ,p) = e- s ^^f^x,b 2 ^ 2 ^ b )b\(z,b 2 ^ 2 ,^)^ (17) 



with 



and 



with 



S0,Q) = - hl [<p\K(b,n b )-J a ^[ lF (g^)-A)~\\n{^ lK {g^))l (18) 



fi{x,b 2 -X F ^)b\{z,b 2 ] CD^i)=e- s ^ Q ' Q ^ft{x,b 2 -QlQ )b\{z,b 2 -QlQ Q ), (19) 



S(b,Q,Q )=-]xi(^\ (K(b,m,) + p ^V(<?(m))) " / 2 ^[7F(5(M);l)-^ln^) 7 K(^))], (20) 
where we have used that 7^ = 7^ to the order in a s considered here. 



1 The author is grateful to John Collins and Ted Rogers for pointing this out. 



B. Pert ur bat ive Sudakov factor 



The various quantities in the Sudakov factor to order a s are given by [10] : 



K(b,/i) = -a 8 ( M ) — H^/Ci) + 0(ai), 

IT 



(21) 
7jc(s00) = 2a s (»)^ + 0(ai), (22) 

7T 

7f(9W, C/A^ 2 ) = a s (^ (I ~ ^ (C/m 2 )) + 0(a 2 s ). (23) 



Here it should be noted that for a running a s , the choice of scale /i, and hence of the integration range over // 
in the Sudakov factor, matters much for the size of the errors here generically denoted by 0(a 2 ). Depending on 
the choice of factorized expression, including the choice of using /j,f, or Q$, the error in the final result for the 
asymmetry may vary considerably in size. We emphasize that the fixed scale Qo choice is not necessarily the 
optimal choice in all cases. 

From these perturbative expressions one obtains the following perturbative Sudakov factors: 
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Including the one-loop running of a s one can perform the /i integrals explicitly. We define 
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such that (dropping non-logarithmic finite terms) 
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(24) 
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(27) 



The above expressions for the Sudakov factor are valid in the perturbative region b < 1/Qo- Strictly speaking, 
at very small b (< 1/Q) the perturbative expressions do not have the correct behavior S — ¥ in the limit 
b — > 0. If this region gives important contributions, it requires modifications of the expressions (for example the 
regularization discussed in [TH]) or else it can lead to artifacts, such as 'Sudakov enhancement' of the cross section. 
For the Sivers asymmetry calculation we find that such a modification is needed if one uses Eq. (15) instead of 



Eq. (16). Using Eq. (16), and hence Eq. (27), leads to only minor contributions from the region b < l/Q, so we 
will use them without modifications. 



C. Nonperturbative Sudakov factor 



As said, the above expressions for the Sudakov factor are valid in the perturbative region b < 1/Qo- Since 
at sub-asymptotic Q values, the Fourier transform involves also the nonperturbative region of large b and we 
explicitly focus on the region Q\ <C Q 2 , we will have to deal with b > l/Qo also. This can be done for instance 
via the introduction of a 6-regulator [19 2 : b — >• b* = b/y/l + & 2 /& 2 nax , such that 6* is always smaller than 6 max . 
One then rewrites W(b) as [15] : 



W(b) = W(K)e- SNp(b \ 



(28) 



An alternative method using a deformed contour in the complex b plane has been put forward in Ref. 1201 . 



where for the function W(b*) a perturbative expression can be used. The nonperturbative Sudakov factor Snp 
to be used here is the one from Aybat and Rogers [TO] 



SNp(b,Q,Q ) 
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(29) 



with 0x = 0.201 GeV 2 , g 2 = 0.184 GeV 2 ,g 3 = -0.13, x = 0.009, Q = 1.6 GeV and 6 max = 1.5 GeV . This form 
is chosen such that at low energy (Q = \/2A GeV) it yields a Gaussian with (kj) = 0.38 GeV 2 that resulted 
from fits to SIDIS data [H] and at high energy it matches onto the form fitted simultaneously to Drell-Yan and 
Z boson production data [22]. Like in AR, x = 0.09 will be chosen, resulting in 



SNp(b,Q,Qo) 



0.184 In 



2Qo 



0.332 



(30) 



Note that Snp has a Q 2 dependence that is in accordance with the phenomenological observation that the average 
partonic transverse momentum grows as the energy (Q or y/s) increases (cf. Fig. 12.3 of [13]). The Q-independcnt 
part of Snp can in general be spin dependent, which means that strictly speaking one should allow for a somewhat 
different Snp in the numerator and denominator of the Sivers asymmetry. At lower Q 2 this can become relevant. 
Although quantitatively the results depend considerably on the choice of Snp, in B01/B09 it was found that the 
Q 2 dependence of the ratio is not very sensitive to it. Below we will briefly comment on it further. 



D. TMDs at small b 



Using the above perturbative and nonperturbative Sudakov factors, we end up with the asymmetry expression: 
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(31) 

Assuming the TMDs are slowly varying as a function of 6* , as they are when the TMDs are taken to be (broad) 
Gaussians, like in B01/B09 3 , yields 



A a b(x,z,Q T ) = 
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for some functions fi(x;Qo), Di(z;Qq) (a priori not coinciding with the collinear parton distribution and frag- 
mentation functions at the scale Qo), fir a (x;Qo), and 
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The approach of AR includes the perturbative expansion of the 6* dependence of the TMDs, which yields In Q^b 2 
terms, but also integrals over momentum fractions and mixing between quark and gluon operators. Clearly this 
is the more sophisticated approach, but it also makes it harder to handle. It was not included in the analysis 
of APR which confronts the Sivers asymmetry evolution with HERMES and COMPASS data. Here we will also 
not include it. The expression for A a b(x, z, Qt) in Eq. (32) thus corresponds to the recent TMD factorization 



based approach discussed in APR. The approach followed in B09 included some 6* dependence beyond that of 
the Sudakov factor, arising from the so-called soft factor, see Sec. [Vj It was found that it makes the asymmetry 
fall off somewhat faster with energy. 

It should be noted that the above simplification of dropping the perturbative tails of the TMDs (the 6* de- 
pendence of the TMDs), every asymmetry involving one fc^-odd TMD will be of the form in Eq. (32), and hence 



proportional to A(Qt)- This has the advantage that the results obtained below also apply to for instance the 
Collins asymmetry in SIDIS. Of course, the latter will be multiplied by a different x and z dependent prefactor. 



3 Note that in B01/B09 the Gaussian width of the Sivers TMD appears in the asymmetry expressions, because of the derivative in 



IV. NUMERICAL STUDY OF THE Q T AND Q DEPENDENCE OF A(Qi 



The expression studied numerically is A(Qt) as given in Eq. ( 
(27) and the nonperturbative Sudakov factor Snp from AR/A 



33), using the perturbative Sudakov factor in Eq. 
R in Eq. Q. In Fig. [I] (left) A{Q T ) is plotted 



for various energies. As can be seen, the asymmetry has a single peak structure, whose magnitude falls off with 
energy. In Fig. fl] (right) A(Qt) evaluated at Qt = Qr.max at which the asymmetry reaches its peak, is plotted as 
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FIG. 1: Left: the TMD factorization based asymmetry factor A(Qt) (in units of M 2 ) at Q - 3.33, 10, 30, 60, 90 GeV. 
Right: the asymmetry factor A(Qt) evaluated at Q-r.max (solid line) plotted as function of Q and compared to a line 
(dotted) with 1/Q ' 68 fall-off, constructed to coincide at the point Q = 30 GeV. 

a function of Q and compared to a very simple power law approximation. This shows that the peak of A(Qt) has 
to good approximation a 1/Q 068 fall-off, which is only slightly faster than the results of B01, where a l/Q 05-06 
dependence was found. 

The reason for focusing on the Q dependence of the peak rather than of the asymmetry at a fixed transverse 
momentum or of an integral of the asymmetry, is simply that for future experiments at higher scales one is first 
of all interested to know the minimal sensitivity required to observe a nonzero asymmetry signal. It is thus most 
interesting to study the evolution in the Qt region where the asymmetry is largest, which is a region that shifts 
towards higher transverse momentum values as the energy increases. It matters less for experimental studies that 
on the sides of the peak where the asymmetry is considerably smaller, it falls off even faster. In Fig. [2] it is shown 
how QT,max nioves towards higher values as Q increases. If this growth of QT,max is not confirmed experimentally, 




100 



FIG. 2: The peak position Qr.max as function Q. 



it likely means that the 6* treatment of the TMDs is oversimplified, i.e. that the perturbative tails of the TMDs 
matter. Also the choice o f Sn p may make a difference. To give an idea of the dependence of the results on Snp' 
multiplying S NP in Eq. Q by a' factor of 2 yields 1/Q 064 and by a factor of 1/2 yields 1/Q 075 . Although 
the power of the fall-off is not affected much, the peak and its position do change considerably, in general by 
40-50% (very similar to what was found in B01). It should be kept in mind though that the S^p used here is 



fitted to available unpolarized data over the entire considered Q 2 range and is therefore very appropriate for the 
denominator of the asymmetry. Variations by a factor of 2 are thus not realistic (remember that they appear in 
an exponent). But what is very well possible is that the Q 2 independent part of Snp in the numerator is smaller 
than what is used here (larger is not allowed by the positivity bound at Q = 2Q ). If it is reduced by a factor 
1/2 it yields 1/Q 079 . Taking into account the uncertainty from Snp, we thus estimate the power to be in the 
range 0.6-0.8. 

V. COMPARISON TO A CS FACTORIZATION BASED APPROACH 

We already compared the results obtained within the recent TMD factorization approach with some results of 
the study B01. Since the approach of B09 was an improved version of B01, also based on the original Collins-Soper 
(CS) factorization, it may be of interest to compare the above results to those that would follow from B09. This 
can be useful for estimating the size of the expected modifications of old results for other asymmetries as well. 

The approach in B09 was based on a fixed scale Qo version of the CS factorization, obtained from Ref. [7j by 
using the relevant renormalization group equations. This resulted in a perturbative Sudakov factor |14j : 
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where the renormalization-scheme-dependent finite term F% = 27^ — 3/2 — In 4 will be dropped. We emphasize that 
although this expression is here denoted by Sbo9, because it was written in this way in B09, it can straightforwardly 
be obtained from the original CS paper by replacing its \i\ — C\jb — » Qo and carrying through all the relevant 
subsequent replacements in the various renormalization group equations. 
This Sudakov factor can be rewritten as 
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which differs from S p (b, Q, Qo) in Eq. (27 1 only by sub-leading logarithmic terms. In the leading (double) logarith- 



mic approximation of fixed coupling constant, all perturbative Sudakov expressions coincide with the well-known 
result of [18. Any numerical difference between the approaches thus indicates sensitivity to single logarithms, like 
from the running of a s . Needless to say, this becomes more important the larger the Q range is in the comparison. 
In the expressions of B09 also a soft term U(b; Qo, a s (Qo)) needs to be included [13] : 



U{b-Qo,a s {Qo)) = l-^Ml CF (\nQ 2 b 2 
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where the finite term F 2 = ln7r + 7^ will be dropped again. All this amounts to inserting in Eq. (31 1: 



f^ a (x,bi;QiQo)D a 1 (z,bi;QiQo) = f& a (x;Q )D<[(z;Qo)U(b*;Qo,a s (Q )), 
/?(*, bl;Q 2 , Qo)D\{z, bl- Q 2 , Q ) - f^x; Q )D b 1 (z; Qo)U(K;Qo, a.(Q )), 



(37) 
(38) 



and replacing S p by 5^09- The approach followed in B09 thus includes some 6* dependence beyond that of the 
Sudakov factor. As said, this leads to a somewhat faster fall-off of the asymmetry with energy. 
Putting this together, the expression to be evaluated is: 



Abo9(Qt) = M 



J dbb 2 7i(hQ T ) #(&»; Qo.MQo)) exp(-S B09 (K,Q,Qo) - S NP (b,Q/Q )) 
f dbbJ (bQ T ) U(b*;Q ,a s (Qo)) exp (-S B0 9(K,Q, Qo) - S NP (b, Q/Q )) ' 



(39) 



In Fig. [3] (left) Abo9(Qt) is plotted for various energies. In Fig. [3] (right) ^4_B09(QT,max) is plotted as a function 
of Q and compared to a simple power law approximation 1/Q 065 (without soft factor the calculation yields a slower 
fall-off: l/Q 058 ). Note that the asymmetry is generally smaller than A(Qt) and that the peak of the asymmetry 
moves more slowly towards higher Qt values. Nevertheless, the power of the fall-off is very comparable to the 
result of the previous section. We therefore expect the estimates of the evolution of the cos 2<fi Collins asymmetry 
discussed in B09 (approximately 1/Q) to not change much when using the recent TMD factorization approach. 
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FIG. 3: Left: the Collins-Soper factorization based asymmetry factor ,4bo9(Qt) (in units of M 2 ) at Q = 
3.33, 10, 30, 60, 90 GeV. Right: the asymmetry factor Aboq (Qt) evaluated at the QT,max plotted as function of Q (solid 
line) and compared to a line (dotted) with a 1/Q 0,58 fall-off, constructed to coincide at the point Q — 30 GeV. 



We note that like previously considered perturbative Sudakov factors, Sbo9 a l so does not have the correct 6^0 
limit. Due to its InQ^b 2 dependence, compared to the InlnQ^b 2 dependence of S p in Eq. (27), the region of very 
small b (below 1/Q) contributes more than before. But it turns out to only matter in the denominator of the 
asymmetry to a modest extent (around QT.max this region contributes 10-20% of the denominator of Abq9(Qt): 
compared to around 5% in the case of A(Qt))- It may be the reason for a somewhat smaller asymmetry. The 
perturbative tails of the TMDs may matter more in this case. These tails contain logarithms of 6*Qo, which only 
become large in the very small b region. Only if the remainder of the integrand has little or no support there, 
these logarithms are of no importance. Given the modest contribution of the very small b region even in this B09 
approach, we will not investigate this issue further here, but it is an argument in favor of using the new TMD 
factorization approach as in A(Qt)- 



VI. CONCLUSIONS 



In this paper the energy scale dependence of the Sivers asymmetry in SIDIS has been investigated within the 
framework of TMD factorization. The perturbatively calculable part of the Sudakov factor is considered at the 
one-loop level, including higher order effects due to the running of the coupling constant, in order to avoid the 
appearance of large logarithms. This study is very similar to the recent one in Ref. [12 (APR), which focussed 
on the low energy region Q < 10 GeV. Here we study a larger Q range and specifically the Q behavior of 
the peak of the Sivers asymmetry, which allows for a direct comparison with earlier results based on Collins- 
Soper factorization in various approximations. Although these treatments differ only beyond the double leading 
logarithmic approximation, the numerical results show that subleading logarithms do matter in the studied Q 
range. The recent TMD factorization based approach indicates that the peak of the Sivers asymmetry falls off 
with Q approximately as 1/Q 07 , somewhat faster than was found in [T3] (1/Q ' 5-0 ' 6 ), but similar to what follows 
from the approach discussed in Q3] which still includes a soft factor (1/Q - 65 ). Since these numerical results 
involve some approximations, such as neglecting the b dependence of the TMDs in the small b region and using 
the same S^p in both numerator and denominator of the asymmetry, the actual fall-off to be determined by 
experiment could be somewhat faster even. From varying the nonperturbative Sudakov factor (also separately in 
the numerator of the asymmetry) we expect a power somewhere in the range 0.6-0.8. Similar moderate differences 
between the approaches based on the first TMD factorization of Collins and Soper (1981) and on the recent TMD 
factorization by Collins (2011) are expected also for other azimuthal asymmetries, such as the Collins effect 
asymmetries studied in [T31 [TJ] . Of course, this conclusion applies specifically to the studied kinematic range. 

In the numerical results the peak of the asymmetry moves towards higher transverse momentum values as 
the energy increases, by approximately 70% over the studied Q range of 3 — 100 GeV. The peak is located at 
a transverse momentum value that is comparable to Qo, where one expects the dominant contribution to the 
b integral to come from the region b ~ 1/Qcb which is the boundary of the perturbative region. Inclusion of 
the perturbative tails of the TMDs may thus affect the location of the peak of the asymmetry. Apart from the 
power of the fall-off in Q, the behavior of the peak as a function of Qt will therefore allow to test the underlying 
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assumptions considered in [12| and in this paper. Clearly, higher Q data on the Sivers asymmetry in SIDIS is 
required to test the evolution resulting from the TMD formalism, which awaits an EIC. 

Note added: upon completion of this paper, a paper on the energy evolution of the Sivers asymmetries appeared 
that addresses related topics [24], in particular questioning the Snp used in APR. As we pointed out, our results 
for the power of the fall off with energy are not very sensitive to the particular S^p used. In addition, a new 
study of the effect of TMD evolution of the Sivers function in semi-inclusive J/ip production appeared recently 
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